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We derive and validate a set of coupled Bloch wave equations for analyzing the reflection and trans¬ 
mission properties of active semiconductor photonic crystal waveguides. In such devices, slow-light 
propagation can be used to enhance the material gain per unit length, enabling, for example, the re¬ 
alization of short optical amplifiers compatible with photonic integration. The coupled wave analysis 
is compared to numerical approaches based on the Fourier modal method and a frequency domain 
finite element technique. The presence of material gain leads to the build-up of a backscattered 
field, which is interpreted as distributed feedback effects or reflection at passive-active interfaces, 
depending on the approach taken. For very large material gain values, the band structure of the 
waveguide is perturbed, and deviations from the simple coupled Bloch wave model are found. 
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I. I. INTRODUCTION 

Photonic crystal (PhC) structures have been proposed 
as a waveguide infrastructure for high-density photonic 
integrated circuits m- Optical amplification is one of 
the fundamental functionalities, required for compensat¬ 
ing attenuation and coupling losses and thus scaling-up 
the number of integrated devices [4j. A major advantage 
of using photonic crystal waveguides, as opposed to stan¬ 
dard ridge-type waveguides, when realizing active struc¬ 
tures with embedded semiconductor gain layers, is the 
use of slow-light effects to enhance light-matter interac¬ 
tions. Thus, the spatial gain coefficient may be increased 
by exploiting slow-light propagation, enabling the real¬ 
ization of short devices suitable for photonic integration 
W- Similarly, nonlinear effects induced by carrier de¬ 
pletion may also be enhanced, leading to devices with 
ultra-low saturation power, which may be of interest for 
optical signal processing applications 0. 

In this work, we theoretically investigate the 
frequency-domain optical propagation properties of an 
active PhC waveguide of finite length embedded in an 
ideal passive periodic PhC waveguide platform as shown 
in Fig. [T] Such structures with site-controlled active gain 
sections may be fabricated using different techniques 0 - 

EE 

As shown in Fig. |T| we approximate the device as a 
photonic crystal heterostucture m with a slowly-varying 
envelope of the imaginary part of the refractive index 
(bottom part), and a fast periodic variation that natu¬ 
rally arises from the penetration of the air holes through 
the gain layer (top part). Such an envelope approxima¬ 
tion in photonic crystal devices is analogous to the treat¬ 
ment of semiconductor opto-electronic devices. Based on 
a perturbative approach [141, effective one-dimensional 
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(ID) coupled-wave analysis has been widely used to in¬ 
vestigate the impact of slow light effects on optical prop¬ 
erties of passive PhC waveguides, e.g. for efficient taper 
design EE Kerr nonlinearities m ee and disorder- 
induced scattering ®EE 



FIG. 1 Illustration of defect photonic crystal waveguide. 
The red part of the structure is active, i.e., this part of the 
membrane structure contains embedded layers of quantum 
wells or quantum dots that may provide gain upon optical 
or electrical pumping. The lattice constant is a and the 
length of the active region is L. 

For perfectly periodic PhC waveguides, neglecting Kerr 
nonlinearity and disorder, a rigorous set of equations for 
the amplitudes of forward and backward propagating un¬ 
perturbed Bloch waves may be derived HOI. In such a 
formulation, the presence of active material in a finite 
section of the PhC waveguide leads to multiple scatter¬ 
ing, which represents material-gain-induced coupling be¬ 
tween the forward and backward Bloch waves of passive 
structures. Such distributed feedback (DFB) effects have 
so far not been considered in slow-light enhanced semi¬ 
conductor optical amplifiers [3 CHH2E We show here 
that it is important to consider the impact of such feed¬ 
back effects when calculating the gain of an active PhC 
waveguide. 
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Alternatively, we may treat the structure as multiple 
waveguide sections, namely an active PhC waveguide sec¬ 
tion interfaced with two semi-infinite passive PhC waveg¬ 
uides on both sides, with distinct sets of Bloch modes 
in the passive and active PhC waveguide sections. In 
this picture, the active-passive interfaces induce reflec¬ 
tions and Fabry-Perot effects are important in determin¬ 
ing the strength of the signal transmitted from the input 
to the output waveguide as well as the back-reflected sig¬ 
nal. We solve the optical multiple scattering problem 
in aggregates of generalized thin films using the Fourier 
modal method (FMM) (or Rigorously coupled-wave anal¬ 
ysis (RCWA)) [HHMI, which has been widely imple¬ 
mented for various periodic photonic structures. Using 
these methods, Bloch waves are computed and sorted in 
both passive and active PhC sections [ 13 - 

In the following, we numerically validate that both ap¬ 
proaches are equivalent in analyzing the proposed active 
PhC devices. A similar discussion for active Bragg grat¬ 
ing structures can be found in [25] , We also analyze the 
situation with very large material gain values, in which 
case the band structure of the photonic crystal structure 
is modified. 

We shall here limit our attention to the simplest type of 
photonic crystal waveguide, with one row of missing holes 
(Wl-waveguide). In that case, slow-light propagation is 
realized close to the Brillouin zone edge and the band¬ 
width over which significant slow-down effects are real¬ 
ized is limited. It was shown, though, that by dispersion 
engineering the waveguide, a slow-light region of signifi¬ 
cant bandwidth can be obtained for a frequency range 
displaced from the band-edge [25.. A systematic ap¬ 
proach for optimizing the group index bandwidth prod¬ 
uct was recently demonstrated [30] . For such dispersion 
engineered structures, the large group index region may 
be centered around an inflection point of the dispersion 
curves and it was shown that strong losses (or gain) can 
induce drastic changes of the dispersion curves [31]. 

The paper is organized as follows: In Section II we 
derive the coupled Bloch wave equations. Section III 
presents the numerical results. We first consider a long 
waveguide with relatively weak gain, and compare the 
coupled wave analysis to numerical results obtained us¬ 
ing the FMM. Secondly we consider a shorter waveguide, 
enabling the consideration of larger absolute gain coef¬ 
ficients as well as the use of a finite element numerical 
technique. Finally, the main conclusions are summarized 
in Section IV. 


II. II. EFFECTIVE ONE-DIMENSIONAL 
COUPLED WAVE ANALYSIS: PERTURBATION 
AND DISTRIBUTED FEEDBACK 


We consider a line-defect photonic crystal waveguide, 
where part of the waveguide is active, as illustrated in 
Fig-S The active material may be implemented as lay¬ 
ers of quantum wells or quantum dots embedded in the 



FIG. 2 Illustration of coupled wave analysis of PhC 
waveguide, (a) Schematic of surface used in the derivation of 
coupled mode equations, with arrows indicating outward 
unit normal vectors to the boundary surface, (b) Schematic 
diagram of coupling between the amplitudes of forward and 
backward propagating Bloch waves induced by the active 
material, which is represented by an imaginary perturbation 
of the refractive index. 


middle of the photonic crystal membrane. In order to 
limit the energy consumption, it is preferable only to 
have active material in the core of the waveguide, as 
demonstrated in j5]. Here, however, we shall consider 
the simpler structure, where the membrane is uniformly 
active in the transverse directions ( xy). This is appropri¬ 
ate for the case of optically pumped structures, as con¬ 
sidered, for example, in [7]. We represent the material 
gain as a weak perturbation to the passive periodic PhC 
waveguide in a finite-length section and use the Bloch 
modes of the passive structure to expand the field of the 
structure including gain. In this formulation, multiple 
scattering in the active PhC structure is represented as 
material-gain-induced coupling between the forward and 
backward propagating Bloch waves of the passive struc¬ 
ture, as illustrated in the top part of Fig. 2jb). This is 
similar to the way disorder is analyzed in [18] . 

In the continuous wave (CW) limit, the Maxwell equa¬ 
tions can be rewritten based on the conjugated form of 
the Lorentz reciprocity theorem mm- 

V • (E X HS 1± + Eq ;± X H) = iccP pert • Eq >± . (1) 

Here, E 0 ,± and H 0j ± are the guided Bloch wave solu¬ 
tions of Maxwell equations for the electrical and magnetic 
fields obtained in the absence of electronic polarization 
perturbations P per t, with ± denoting forward/backward 
propagating fields. The fields E and H are the self- 
consistent solutions in the presence of electronic polar¬ 
ization perturbations, which arise due to the presence of 
active material, e.g. leading to stimulated emission and 
thus gain of the propagating fields. We may represent 
the unperturbed Bloch waves as follows: 

Eo,± = ^e±(x,y, 2 ) exp(±ifc 2 z) exp(—iwt), (2) 

H 0 ,± = ^h ± (ar, y, z) exp(±ifc, 2 :) exp(-iwt), (3) 

with e±(x,y,z) and h±(x, y,z) being the normalized 
complex amplitude Bloch functions, ui the frequency and 
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k z the Bloch wave number along the length of the waveg¬ 
uide. 

The unit rms power flux P z and unit rms electric and 
magnetic stored energy W in a supercell are given as 

Pz = \ / [5?{e x h*} • z] dS, (4) 

J s 

w = \J v [eo^ft(r)|e| 2 + Mo|h| 2 ] dH = °^P Z . (5) 

Here, n gz = c/v gz is the group index, with c and 
v gz = duj/dk z being the speed of light and the group 
velocity, £o and /io are the vacuum permittivity and per¬ 
meability, rib(r) is the background refractive index, a is 
the lattice constant of the photonic crystal, S indicates 
the transverse plane at position z, and V is the volume 
of a PhC supercell. 

In the limit of a weak electronic perturbation, we ex¬ 
press the self-consistent solutions as superpositions of the 
guided Bloch waves for the passive structure, with slowly 
varying amplitudes t/j±(z): 

E ~ ip + (z)F, 0 ,+ + tM~) E o+ c.c., (6) 

H~^ + ( 2 )Ho, + + V--(^)H 0 ,_+c.c.. (7) 

By using Eqs. 0, 0, 0 and 0, Eq. (jTJ) can be for¬ 
mulated as a set of continuity equations for the forward 
and backward propagating amplitudes: 

± iv • [V»±(*)K{e X h*}] = iwP pert • E* )± . (8) 

Eq. 0 can be derived using the divergence theorem for 
a closed surface S enclosing the simulation domain: 

ty±(z)St{e x h*} ■ n] dS = iw [ (P pe rfES i± ) dV. 

Z s Jv 

(9) 

Here, n is the outward unit normal vector of the closed 
surface depicted in Fig. 0a). This closed surface can be 
chosen conveniently to comply with the experimental ge¬ 
ometry considered (near field measurement ]33l or trans¬ 
mission measurement ; 53j), as well as the simulation do¬ 
main and boundary conditions (BCs) (e.g. periodic BCs 
[55] or perfectly-matched layer (PML) BCs [25]). The 
left-hand side (LHS) of Eq. 0 describes the possible 
channels for the power flux in/out of the surface and the 
right-hand side (RHS) describes the source term. Eq. 0 
is in the form of a set of implicit integral equations. As¬ 
suming that the considered modes are transversely con¬ 
fined, the surface integral (LHS of Eq. 0) reduces to two 
integrals between the planes z = 0 and z = L, described 
as: 


± \<ff> s [ip±(z)$l{e x h*} • n] dS 1 

~ ± fy d z tp±(z) dzx \ f s [5i{e x h*} • z] d S. (10) 

For simplicity, we describe the carrier-induced complex 
susceptibility perturbation as a product of a complex con¬ 
stant Xpert and an active material distribution function 


F(r): 


Ppert ~ E-OXpert A(l )E, 


( 11 ) 


where F(r) = 1 (= 0) in the active (passive) region. 

We now obtain the coupled-mode equations from 
Eq. 0 with substitution from Eqs. ([T()|)-(|lT]): 


d z ip + (z) = —n gz x P ert [S(z)t/j + + n*(z)e l2kzZ ipJ\ (12) 


d z tjj_(z) = — l —n gz x P ert [8(z)ip- + n(z)e 


i2 k z z 


with 


S{z) = 


k(z) = 


a fs [ g o F ( r )l e l 2 ] dS 
fv [eonb(r)|e| 2 + ^o|h| 2 ] dV’ 
a f s [£ 0 E(r)e • e] d S 


V’+](13) 


(14) 


(15) 


fv [ e on 2 (r)|e| 2 + Po|h| 2 ] dW 

Here, S(z) and k(z) denote the complex-valued propaga¬ 
tion and backscattering coefficient induced by the pertur¬ 
bation due to the active material. The above ID coupled 
wave equations are mathematically equivalent to a 2 x 2 
scattering matrix problem in a stack of thin films (at each 
z coordinate) and readily solved numerically. 

We note the similarity of the above coupled mode equa¬ 
tions with the coupled mode equations originally derived 
by Kogelnik and Shank for distributed feedback lasers 
[5B1 . However, while those equations describe the slow 
variation of the complex amplitudes of plane waves (see 
also m for application to PhC waveguides), Eqs. (H0- 
(13) describe the slowly varying amplitudes of the Bloch 


waves of the corresponding passive PhC. This has the 
consequence that in the absence of a polarization pertur¬ 
bation, i.e., Xpert = 0> the forward and backward Bloch 
waves do not couple but remain uni-directional propagat¬ 
ing waves, as required. The presence of a perturbation, 
however, leads to coupling of the Bloch waves, since the 
material perturbation gives rise to distributed feedback. 
This is analogous to the scattering between Bloch waves 
induced by structural disorder [18] , although in that case 
the backscattering sites are randomly distributed. 

In general, the harmonic term exp(±i2fc,2:) leads to 
relatively low backscattering efficiency due to phase mis¬ 
match between forward and backward propagating Bloch 
waves. Propagation losses decrease the effective propaga¬ 
tion length and further damp the backscattering. In pas¬ 
sive PhC waveguides, one can integrate similar coupled¬ 
mode equations approximately by neglecting the dis¬ 
tributed feedback terms as long as the backscattered 
power is weak m On the other hand, in the case of 
positive gain, the effective length, over which wave inter¬ 
action takes place, increases, and a significant backscat¬ 
tered wave may be built up, as we shall show here. 

Eq. (12)-(13) must be supplemented by boundary con¬ 
ditions at the passive-active interfaces. The boundary 
conditions [38: are: 

ip+{0)=r 1 'tp-(0) + ipo, (16) 

4>-(L)exp(-ik z L) = r 2 ip+(L) exp(ik z L), (17) 
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where rq and r 2 are the amplitude reflectivities of the left 
and right passive-active interfaces and ipo = 1 indicates 
unity incident field from the left interface. 

Finally, we note that the coupled wave approach pre¬ 
sented here permits the inclusion of disorder-induced 
losses pSJ US] and carrier dynamics [8j, which is beyond 
the scope of this paper. 


III. III. SIMULATION RESULTS 

For the numerical investigations, we consider a two- 
dimensional (2D) triangular lattice of air holes (hole ra¬ 
dius r = 0.25a) embedded in a semiconductor membrane 
with a dielectric constant of n% = 12.1. The waveg¬ 
uide is a W1 line-defect with a single row of air holes 
omitted. The presence of material gain or absorption 
is modeled via an imaginary refractive index perturba¬ 
tion rii, corresponding to a material gain g 0 = — 2rijw/c. 
The corresponding susceptibility perturbation is Xpert = 
—nj + i2nbrii. In practice, the value of the gain coefficient 
is controlled via the charge carrier density. We notice 
that the ratio between the real and imaginary part of the 
susceptibility change is close to zero, corresponding to a 
near-zero value for the linewidth enhancement factor. In 
the results to be presented, we shall quantify the gain via 
the imaginary part of the refractive index; as examples, 
rii = —10” 3 corresponds to a material gain of go = 81.1 
cm” 1 at Ao = 1550 nm and go = 96.7 cm” 1 at Ao = 1300 
nm. 

In the following we separately consider two waveguide 
lengths in order to illustrate and discuss two different 
physical effects. Section A considers the case of rela¬ 
tively long waveguides, but weak material gain, where a 
strong backscattered signal may build up. This situation 
will be shown to be accurately modeled by the derived 
coupled Bloch mode equations. Section B considers short 
waveguides, but with strong material gain, in which case 
the perturbative treatment breaks down and the imag¬ 
inary contribution to the refractive index changes the 
band structure. 


A. A. PhC waveguides with a long active section 

The length of the active PhC waveguide considered in 
this section is L = 100a. Numerically, such long waveg¬ 
uides are conveniently treated using the FMM that ex¬ 
ploits the periodicity for obtaining the Bloch modes in 
both passive and active PhC sections. Specifically, it 
suffices to analyze a single unit cell (of length a), the 
eigensolutions of which are the Bloch modes [25, 27], 
and the propagation through the full length (L = 100a) 
is afterwards handled analytically via Bloch’s theorem. 
Active-passive reflection and transmission is handled us¬ 
ing a scattering matrix approach [39]. 

Figure [3] shows the transmitted and reflected power 
flux of an input CW light exciting a slow light mode 



Imaginary refractive index, -n. 


FIG. 3 Power reflection and transmission for the PhC 
structure in Fig. |l]as functions of the imaginary part of the 
refractive index for a slow light mode at oja/(2nc) = 0.2064 
with a group index n gz ~ 25 (see Figs. |4](a) and (b)). The 
length of the active waveguide is L = 100a. Solid lines are 
results from the ID coupled wave analysis, with 
|n| 2 = |t- 2 | 2 = 0. Markers are 2D FMM simulation results. 


(wa/( 27 tc) = 0.2064, n g ~ 25, see Figs. E]]a) and (b)) as 
functions of the imaginary part of the refractive index. 

The markers show the results based on 2D FMM sim¬ 
ulations, while the solid lines are the results based on the 
solution of the ID coupled Bloch wave equations. As is 
apparent, the two approaches agree well. In both cases, 
a non-zero value of the imaginary part of the refractive 
index leads to the build-up of a reflected signal. The 
dashed line shows the ideal slow-light enhanced device 
gain/absorption [7], G = exp(2n gz ribgo f Q S(z)dz) = 
exp(T^-g 0 L), r being the confinement factor, and the 
slope of which is proportional to the group index. It is 
seen that the simple slow-light enhanced gain expression 
in this case well accounts for the numerical results for 
absolute gain values up to about 25 dB. Beyond that the 
transmission starts to decrease, due to the build-up of a 
strong backreflected signal. 

In the absorption regime (—71* < 0), the increase of the 
material absorption leads to an increase of backscattered 
power flux, but the weak backscattered signal damped 
by material loss has negligible effect on the transmit¬ 
ted power. On the other hand, in the regime of positive 
gain (—m > 0), the forward as well as backscattered 
signals are amplified and the effective length over which 
the waves interact is increased. As the material gain in¬ 
creases, the transmitted power reaches a maximum value 
and then steeply decreases. 

We now briefly discuss the interpretation of the results 
in the context of the two different numerical techniques 
used, which effectively correspond to the use of different 
basis sets for representing the full solution. In the cou- 
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FIG. 4 Propagation characteristics of a line defect photonic crystal waveguide with a gain section of finite length, L = 100a. 
(a) Band diagram and (b) group index of passive PhC waveguide Bloch mode, (c) Power transmission and (d) reflection as 
function of the imaginary part of the refractive index and frequency. The results in (c) and (d) are obtained by numerically 
solving the ID coupled Bloch wave equations. 


pled Bloch mode approach, the total field is expanded on 
the Bloch modes of the passive waveguide (see Eqs. i> 
and §)• The inclusion of material gain or absorption 
corresponds to a periodic perturbation of the imaginary 
part of the refractive index and induces coupling between 
the passive Bloch modes (DFB) that are no longer solu¬ 
tions to the full, perturbed problem. In the FMM ap¬ 
proach, the local solutions are instead expanded on the 
active Bloch modes, i.e., the solutions when including the 
imaginary index perturbation. For an infinitely long ac¬ 
tive waveguide, these solutions are uni-directional Bloch 
modes with exponentially growing powers, which renders 
this concept unphysical. In the case of a finite-length ac¬ 
tive waveguide, the appearance of a backreflected signal 
occurs due to the mismatch of the Bloch modes of the 
passive and active parts of the waveguide. This mode 
mismatch leads to finite interface reflection coefficient 
and the appearance of Fabry-Perot (FP) resonances that 
depend on the active region length. The equivalence of 
the DFB and FP pictures was previously discussed in the 
context of simple ID optical grating structures [25] , 

Figures [4|c) and (d) show the transmission and re¬ 
flection as a function of frequency for varying imaginary 
part of the refractive index. For a given frequency, as the 
imaginary part of the refractive index (material gain) in¬ 
creases, both transmission and reflection increase, reach 
a maximum and then start to decrease. For certain com¬ 
binations of frequency and material gain both the trans¬ 
mitted and reflected beam diverge. These points corre¬ 
spond to the fulfillment of the lasing oscillation condition 
for a laser with length given by the active region. The 
corresponding frequencies correspond to discrete values 
of the wavenumber, with the real part of the wavenumber 


fulfilling the condition 

^-Re(fcz)^ 2 L = 2'Km, m = integer. (18) 

Taking into account saturation of the gain medium due 
to stimulated emission, the divergences correspond to the 
onset of self-sustained lasing, with the imaginary part of 
the refractive index being clamped to its threshold value. 
Again, we notice the similarity of these results with the 
classical analyses of distributed feedback lasers (3B; and 
active Bragg gratings [28] . 

For frequencies closer to the bandgap, the resonances 
appear at smaller gain values, consistent with the en¬ 
hancement of the gain per unit length due to slow light 
propagation [7|. Notice that the gain per unit time in¬ 
side the material is not changed, with the consequence, 
for lasers, that it is only the contribution of mirror losses 
to the laser threshold that is decreased when exploiting 
a slow-light mode. 

B. B. PhC waveguides with a short active section 

In this section we consider a short active PhC waveg¬ 
uide section, of length L = 20a, but allow for higher val¬ 
ues of the material gain coefficient in order to explore the 
regime investigated in [23]. As the size of 2D simulation 
domain decreases, it becomes feasible to solve the peri¬ 
odic optical waveguiding problem [40] in the frequency- 
domain using the finite element method (FEM) with nu¬ 
merically exact Bloch mode excitation and/or absorptive 
boundary conditions [41) . 

Figure [5] shows the transmitted and reflected power 
flux of an input CW light exciting a slow light mode 
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FIG. 5 Power transmission and reflection in (a) dB scale and (b) linear scale as functions of the imaginary part of the 
refractive index, for a short active waveguide of length L = 20a. The frequency is taa/(27rc) = 0.2075, with a corresponding 
group index of n gz ~ 18.7. Solid and dash-dotted lines are ID coupled wave analysis results with zero passive-active interface 
reflections (n = r 2 = 0). Markers are 2D FEM simulation results. 


(wa/(27rc) = 0.2075, n gz ~ 18.7) as a function of the 
imaginary part of the refractive index. The solid and 
dash-dotted lines are obtained from the ID coupled wave 
equations, while the markers are corresponding numeri¬ 
cal results based on 2D FEM. The results are shown in 
both logarithmic and linear scales in order to appreciate 
the differences between the two approaches. We see that 
the coupled wave analysis reproduces most features seen 
in the full numerical FEM simulations for values of ma¬ 
terial gain/absorption coefficient up to 4000 cm -1 (\m\ 
up to 0.05) at a wavelength of 1550 ran, which is large 
compared to gain values typically obtained in semicon¬ 
ductors. For smaller gain values, the shorter device shows 
characteristics similar to the longer device, with the abso¬ 
lute gain being smaller for the same gain coefficient. We 
attribute the difference between the coupled wave model 
and the full numerical result to a change of the effective 
gain coefficient appearing as the pre-factor of T + and 
4U in the RHS of Eqs. |l2|-(fT3|). Thus, for a large ab¬ 


solute value of the imaginary part of the refractive index 
(large gain or absorption), the band structure is modi¬ 
fied such that the group index is reduced and thereby 
also the slow-light enhanced gain and absorption. This 
is in agreement with the result obtained in [23] ■ where 
the band structure of an infinitely periodic structure was 
analyzed. In that case, there is no effect of backscat- 
tering, since a uni-directional Bloch wave of the active 
structure is considered, but consistent with the results 
presented here the slow-light enhancement was observed 
to decrease with the absolute value of the material gain 
coefficient. 


We interpret the gain-induced change of effective group 
index (slow-light enhancement factor) as follows. The 
origin of the photonic crystal bandgap is the destruc¬ 
tive interference of multiply scattered waves. However, 


in the presence of gain (or absorption), the destructive 
interference is incomplete. For a defect photonic crys¬ 
tal waveguide, this means that the band-edge is smeared 
and the group index no longer diverges when approaching 
the band-edge. In consequence, there is a gain-induced 
reduction of the group index, and thereby a gain-induced 
reduction of the slow-light enhancement of light-matter 
interaction. 

In the example above we considered a relatively modest 
group index of ~ 18.7, but when approaching the band- 
edge further, the (passive) group index rapidly increases, 
and the saturation effect discussed above becomes more 
apparent, i.e., sets in at lower values of the absolute gain 
coefficient. 


IV. IV. CONCLUSIONS 

We have derived a set of coupled Bloch wave equa¬ 
tions that are useful for simulating and understanding 
the properties of active photonic crystal waveguides. The 
coupled Bloch wave model explains how slow-light effects 
can be used to enhance the absolute gain of a photonic 
crystal waveguide of given length, an effect recently ob¬ 
served experimentally j7]. The model also shows, how¬ 
ever, that the presence of material gain (an imaginary re¬ 
fractive index perturbation) induces distributed feedback 
effects, additional to the slow light enhancement of modal 
gain, which leads to the build-up of a backscattered field 
that limits the attainable gain. The coupled Bloch wave 
model is compared to numerical simulations based on 
the Fourier modal method as well as a frequency domain 
finite element solver. In general, very good agreement 
between the different numerical approaches is found, al¬ 
though, depending on the approach used, the presence 
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of a backscattered field calls for different, but consistent, 
interpretations. For very large values of the material gain 
(or absorption) coefficient, the band structure of the ac¬ 
tive photonic crystal itself is modified, leading to further 
changes of the slow-light enhancement, which is an effect 
not captured by the ID coupled Bloch wave analysis. 

Finally, we notice that the transparent and computa¬ 
tionally efficient formulation of the coupled Bloch wave 
equations presented here allows for the inclusion of ad¬ 
ditional relevant effects. In [18], random disorder effects 
in passive structures were analyzed using a similar cou¬ 
pled Bloch wave analysis. With the recent observation of 


lasing in Anderson-localized modes due to disorder in ac¬ 
tive photonic crystal waveguides [42], it appears natural 
to combine the two approaches. 
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